Effective fiber-coupling of entangled photons for quantum communication 
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We report on theoretical and experimental demonstration of high-efficiency 
coupling of two-photon entangled states produced in the nonlinear process of 
spontaneous parametric down conversion into a single-mode fiber. We deter- 
mine constraints for the optimal coupling parameters. This result is crucial for 
practical implementation of quantum key distribution protocols with entangled 
states. 
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Entangled-photon pairs generated in the nonlinear pro- 
cess of spontaneous parametric down conversion (SPDC) 
are proven to be a highly desirable meansi for practical 
quantum cryptography^. The main difficulty of practical 
utilization of such system usually stems from a relatively 
low photon collection efficiency because of the complex 
spatial distribution of SPDC radiation and due to the 
broad spectral width of entangled-photon wave packets. 

The problem of coupling entangled photons into a fiber 
has been considered before by Kurtsiefer et al, 3 . As- 
suming the pump to be a plane wave the emission angle 
of the SPDC has been calculated as a function of the 
wavelength. The waist of the focused pump beam has 
been chosen to maximally overlap the "impression" of 
the Gaussian mode of a single-mode fiber on the crys- 
tal^. It has been pointed out that the coupling efficiency 
may be significantly affected by transverse walk-off. 

In this letter we present a significantly modified ap- 
proach allowing us to achieve a high-efficiency coupling 
of the SPDC pairs into single mode fibers. In particu- 
lar, we demonstrate how the pump beam waist, crystal 
length, optical system magnification and the fiber mode 
field diameter (MFD) must obey a precise joint relation 
in order to ensure high coupling efficiency. We describe a 
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specific model that allows us to determine a scaling law 
accounting for all the real experimental parameters. We 
demonstrate that obtained experimental data are in good 
agreement with the proposed model. 

The two-photons state can be written as^ 

|*>SPDC = J duj dq D J duj e dq e $(q , uj ; q e , oj e ) ■ 

•aj;(w ,q )d];(w ,q e ) |0). (1) 

The function $(q , uj ; q e , u e ) = £ p (q D + q e , uj + 
w e )x^ 2 H ( lo,w ;q e ,u; e ) accounts for the phase match- 
ing conditions. £p(') represents the amplitude of 
the plane-wave expansion of the pump field and 
X (2) (qo,w ;q e ,w e ) = j° L dz exp[-i Afc(q D , w G ; q e , w e ) z] 
with Ak = k p — k oz — k ez . Inside the crystal the z- 
component of the wave- vector is defined as k z (q, u) — 
y/[u>n(q, <jj)/c] 2 — |q| 2 . All the information on the state 
is given by the amplitude Ai,2(xi,ii;X2,£2)^ of de- 
tecting the SPDC two-photons in space-time events 
at (xi,t\) and (X2,i2)- The Fourier transform with 
respect to t\ and t% of the two-photon amplitude 
-4i,2(xi,£i;x2,f 2 ) is given by 

-4i,2(xi,w ;x 2l £4; e ) = / ciq G dq e $(q Q , lo ; q e , u e ) ■ 

•7ti(xi;q ,w )?i2(x2;q e ,We) (2) 



TLj (xj ; q, u>) (j = 1,2) being the Fourier transform of im- 
pulse response functions hj(xj,x;uj) of the optical sys- 
tems^ through which the two photons propagate from 
the output face of the crystal to the detection plane. 

The coupling of the photon pairs into fibers can be 
considered a problem of maximizing the overlap between 
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the two-photon amplitude A1.2 of entangled-photon state 
in the detector plane with the field mode profiles of 
single-mode fibers^. Assuming a quasi-monochromatic 
and quasi-plane wave travelling in the z-direction and 
defined on the two-dimensional continuous space of the 
detector planes, we can express the electromagnetic field 
operator 2?(+) (x, u>) in terms of a linear superposition of 
electromagnetic field operators, ci m (w) and Ck{oS) asso- 
ciated with a complete orthonormal set of functions^. 
Choosing conveniently the noncontinuous set of guided 
modes ipi m (x,uj) and the continuous set of the radiation 
modes ?/>fe(x, cj), respectively^, the field operator can be 
decomposed as l' + '(x,ii;) = J2i, m ^m(x, w) ci m (ui) + 
J dkipk(x.,uj) Ck(uj) where the new operators c a {oj) are 
defined as c a {uj) = J dx-0a(x, w) E^ + \x, w) with a — 
(7,m) or k, and obey the usual bosonic commutation re- 
lations [c a (ui), Cp(uj)] — 5 a _/3. The amplitude in Eq. @ 
can be expanded in terms of guided and radiation modes. 
The coefficients of the expansion for two guided modes, 
(l,m) and (I , ra ), are given by 
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tmj'm 1 = / dXl dX2 ^l,2(xi, X 2 ,W e ) • 

■^ l 2\^o)^{^ e ) (3) 

Coupling into a single-mode fiber can be quantified by 
the coupling efficiency parameter r]{ c defined as the ra- 
tio of the probability to find two photons in the guided 
modes over the square root of the product of the proba- 
bility to find one photon in a guided mode independently 
of the detection of the other photon2ii£. In the case of 
a single-mode fiber, when only the linearly-polarized fun- 
damental mode LPqi ipQi is allowed^, the coupling effi- 
ciency takes the simple form 
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y/pW . 75(2) 
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The numerator of this expression is given by 
pi 1 '*) = J dbj doje I -4oi',oi( too, uJ e ) | 2 where -4oi,oi ( w ° i w e ) 
is given by Eq. The contributions at the 

denominator of Eq. Q are given by V 1 - 1 ) = 

J duj duj e J dx 2 \J dxiA !2 (xi,w ;x 2 ,w e ) Voi^^i^o)! 2 
and analogous expression for P^ 2 \ The maximum 
coupling, i.e. rjf c — 1, is reached when the two-photon 
amplitude in Eq. (J2J is the product of single-mode field 
profiles of two fibers. 

We now consider a model that includes propagation of 
both fields through two equal infinite ideal lenses with- 
out an aperture limit. We also assume that the output 
plane of the crystal, at a distance dbi from the lenses, is 
imaged on the fiber plane, at a distance d a i from the 
lenses, i.e. l/d bl + l/d a \ — 1/f. The amplitude © 
becomes ^4i j2 (x!, uj ; x 2 , uj e ) oc <&(/z.Xi, uj ; /ix 2 , u> e ) where 
ti = dh\/d&\ = db\/ f — I is the inverse of the magnifica- 
tion and ^(x', ui ; x", uj e ) is the 2-D inverse Fourier trans- 
form of the matching function $(q , lo ] q e , u> e ), which 



we calculate in the paraxial and quasi monochromatic 
approximation. Inside the crystal, such approximations 
allow us to consider only first terms in the exponential 
expansion of expression x^ 2 ^(q ,o; ; q e , ui e ). For a type- 
II non-collinear configuration and assuming a pump field 
to be factorable in terms of frequency and wave vectors 
£ p (q +q e , Up) — £ p (ujp) E p (q +q e ), we can calculate the 
inverse Fourier transform with respect to the frequencies 
and obtain $(x' , T + t/2; x" , T - t/2) ~ Il DL (t)S p (T- 
At/D)E p [(x' +x" - At/D)/2]6(x -x" -Bt/D) where 
we have introduced t = t\ — £ 2 and T = (ti + t2)/2. The 
function HDL(t) has value 1 for < t < DL and zero else- 
where. The vectors A = 2M p - M, B = M + 2Q/K de- 
pends on M and M p , which are the spatial walk-off vec- 
tors for the extraordinary fieldsii at the generation and 
pump frequency, respectively. Q is the transverse wave- 
vector associated with perfect phase matching along the 
intersection of cones, and K — 2K K e /(K + K e ) rep- 
resents a mean value of wave- vector for generated pho- 
tons inside the crystal. We have also introduced D — 
l/u — l/u e and A = l/u p — (l/2u G + l/2u e ) where Uj 
is the group velocity for the j-polarization. The expres- 
sion has a simple physical meaning. Due to the local- 
ity of the interaction, as dictated by the delta function, 
photons are created in pairs at each point of the crystal 
illuminated by the pump field, E p (x). After their birth, 
photons propagate in the dispersive nonlinear crystal en- 
vironment experiencing longitudinal D and spatial walk- 
off, M. They spread with respect to each other in time 
and in transverse direction, according to the travel dis- 
tance, z. The transverse spread contributes through two 
distinct processes when multiplied by the crystal length 
L. The product \A\L represents the shift between the 
generated pairs and the pump field, which is also ex- 
traordinary and hence has a walk-off. |B|i is the spread 
between the two entangled photons generated from the 
same pump photon. This vector contains a contribution 
due to a spatial-walk-off, |M|L, and one more term , 
2|Q| L/K, which represents the transverse distance be- 
tween pairs generated at the input face of the crystal 
with respect to the ones generated at the output face, 
due to the geometry of optical propagation inside the 
crystal. If we assume a pump beam with Gaussian pro- 
file, E p (x) = exp(— |x| 2 /2r 2 ), as well as a mode field 



profile, V'j(x) = exp(— |x| /2w ) /i/ttw, we can derive a 
closed expression for the coupling efficiency, namely 
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(l+e) erf(<r c ) 



(T 2 



erf (0-1 ) erf(cr 2 ) 



where £ = w/j,/r p and 



(5) 



(6) 



The parameters a x = \M p \ 2 + \Q\ 2 /K 2 , a 2 = \M P - 
M| 2 +|Q| 2 /-fT 2 , and [3 = |M| 2 +4|Q| 2 /A" 2 are determined 
only by the crystal parameters and phase matching ge- 
ometry. The expression for r)f c depends on the size of 



3 



fibers effectively imaged onto the crystal plane, wfi, and 
the crystal length, L, scaled to the pump beam waist into 
the crystal, r p . 

The experimental verification of the above model 
has been performed using an actively mode-locked 
Ti:Sapphire laser, which emitted pulses of light at 
830 nm. After a second harmonic generator, a 100-fsec 
pulse (FWHM) was produced at 415 nm, with a repeti- 
tion rate of 76 MHz and an average power of 200 mW. 
The UV-pump radiation was focused to a beam diameter 
of 150 /mi inside a BBO crystal cut for a non-collinear 
type-II phase-matching using a f=75 cm quartz lens (see 
Figure P). 
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Fig. 1. Sketch of the fiber coupling experiment. 

Two points of intersection of the ordinary and ex- 
traordinary conesi^ (3.5 deg relative to the pump direc- 
tion) were imaged into a single mode fiber with a MFD 
of 4.2 /im using two f=15.4 mm coupling lenses placed 
78 cm from the output face of the non-linear crystal. 
We obtain a coupling efficiency of ~ 18% with a crys- 
tal of 3 mm length. We achieved ~ 29% coupling effi- 
ciency using a 1 mm crystal in the same experimental 
setup. Taking into account the approximately 50% de- 
tector efficiency and 85% transmittance of filters, one 
can determine that effective fiber coupling coefficients 
can reach 42% and 68%, respectively (same results have 
been obtained with both 11 nm bandwidth filters cen- 
tered at 830 nm and pass-band filters). This result is 
illustrated in Figure |U where a solid line indicates the- 
oretical curve obtained from our model Eq. J5J. The 
parameters used have been calculated for crystal dis- 
persion at 830 nm: |M p | = 0.07631, |M| = 0.07243, 
and Q/K = |M|/2 = 0.036215. The fiber parameter 
is w — MFD/2V2 - 1.48 /tm, the pump beam radius 
r p ~ 53 /im, and /i ~ 49. 

Examining the dependence of r]f c vs. crystal length 
(see Figure 01 for different values of parameter n one can 
notice, for example, that one cannot reach a coupling 
efficiency greater than ~ 50% for r p ~ 53 /tm and for 
crystals longer than 2 mm. 



Fig. 2. Coupling efficiency rjf c from Eq. J^J vs. the crys- 
tal length L for different values of parameter fi. The 
solid line indicates the theoretical curve obtained using 
the parameters described in the text. The two points 
correspond to the experimental data. 

In conclusion, we have evaluated the dependence of 
coupling efficiency of SPDC pairs into single-mode fiber 
on several major experimental parameters. We obtained 
an analytical expression for the coupling efficiency r/f c 
that allows us to characterize the importance of the spa- 
tial walk-off and to choose appropriate values of experi- 
mental parameters to reach high coupling efficiency. Nu- 
merical inspection of the general expression Q with a 
more sophisticated model of the impulse response func- 
tion of our system can allow us to further increase the 
coupling efficiency. Acknowledgements This work was 
supported by MIUR (Project 67679). G.D.G and A.V.S. 
also acknowledge support by DARPA and NSF. 
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